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We introduce a dissipative Ermakov-Pinney system with the dissipation based on the corresponding 
integrable Abel equation. Based on a particular solution, the general solution to dissipative 
Ermakov-Pinney equations of this type is constructed, including the cases of higher-order Reid 
nonlinearities, by using the Ermakov invariant of the system. 
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I. INTRODUCTION 

The nonlinear Ermakov-Pinney (EP) equation has long been known to have profound connections with the linear 
i-C ' equations of identical operatorial form without the inverse cubic nonlinearity and because of this it has been seen as an 
example of 'nonlinearity from linearity' [l|. A historical overview has been written by Leach and Andriopoulos and 
the fundamental importance of the EP equation for parametric oscillators, both classical and quantum-mechanical, 
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with their vast application reaches, is well established in the literature. 
As known, in 1880, Ermakov introduced the following system 



jj + uj 2 (t)y = , 

X + LU 2 (t)x — KX 

which is endowed with the so-called Ermakov dynamic invariant 



.'• H - j)~(t)x - kx 3 = 0, k>0, ^ 



> 

^sO ' fV\ 2 

lO ' Io-K(x,y,x,y) = n I - + {xy-yx) . (2) 

H : . 

We label here this invariant by the tuning parameters of the nonlinear term of the two equations. 
A century later, Ray and Reid 4] and also Sarlet [|[ introduced generalized Ermakov systems and their associated 
invariants of the type If, g , where / and g are functions of the quotient variable y/x. Recent works on the connections 
t-H ' of the Ermakov systems with the nonlinear superposition principle belong to Carihena and collaborators [6| . When 
/ = g = —k and one knows a particular solution x p , the invariant can be written in the following form 
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In this paper, we start by shortly presenting the factorization of linear dissipative ODE which we use to introduce 
a special EP equation with dissipation determined by an integrability condition for the Abel equation of the first 
kind Q. Next, we include this special Ermakov-Pinney equation in a new type of Ermakov system and determine its 
general solution by means of the attached invariant. 

II. LINEAR DISSIPATIVE ODES AND THEIR FACTORIZATION 

Linear dissipative ODEs of the type 

H< + s(CK + h ((> = o (4) 
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can be factored similarly to their nonlinear counterparts Q, i.e 



- *x ] « = , (5) 



which leads to 



Furthermore, one gets 



where the prime stands for 
Then, we have 



«CC - ^ ( $ i(0 u ) - $ 2"c + $ i $ 2« = . (6) 



HC - (*i + *2)uc + (*i*2 - = , (7) 



3(C)- -(*i + « 2 ) 

HO = $i$ 2 - *i , 

which is a system of equations slightly different from that first used by Rosu and Cornejo-Perez [8[ for nonlinear 
equations of the same structure as Q. This system leads immediately to the following Riccati equation 

$; = _ _ h ( 8 ) 

Moreover, since we can always reduce the dissipative ODE Q to a canonical form without dissipation via transfor- 
mation u(() = u(()e~ 2 /f(C) d C 5 then we will obtain 



cc + h(()u = 0, (9) 



where = h(C) — \g 2 {0 — \ then, without loss of generality, we will always solve the Riccati equation 



with g(C) = 0. 

From ([5]), we also have 

U{ = $iu (10) 

which leads to 

- = 4MC (ii) 

u 

from where 

u = c 1 ef C * 1 W dx . (12) 

In other words, for all functions h(C) for which the Riccati equation is integrable, one first determines the factoring 
function $ 1; and writes the particular solution of (QJ as the exponential in (|12p . 

III. EP EQUATION AS A LINEAR DISSIPATIVE ODE 

Take now the EP equation: 



cc+p(C)« + « =o. ( 13 ) 



and consider it as a linear dissipative ODE of the type (U]), so <?(£) = and h(Q = p(£) + cu 4 . The corresponding 
Riccati equation is 

$[ = -$2 _ h ( 14 ) 
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and from (|T2|) , one can write 

h(0= P (Q+c(c 1 ef C *^ d *y 4 . (15) 

As is well known, if we have two linear independent solutions u\, and ui to 

u a + h(()u = 0, (16) 
then the solution of the corresponding EP equation 

v cc + h(()v + cv~ 3 = (17) 

is given by Q 
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v ® = H-w> (18) 

where W is the Wronskian of the two linearly independent solutions. 

A particular solution to fp~6|) that will be later used in finding a general solution can be obtained as follows: Let's 
take /i(£) = A 2 > 0, a constant, which is a case corresponding to the simple harmonic oscillator equation 

m cc + A 2 u = 0, (19) 

with solutions u\ = cosA£ and U2 — sinA£ and Wronskian W = A. 
The EP equation for this case reads 

v (c + X 2 v + cv~ 3 = (20) 

and its solution can be written in the form 



1 - (1 + sin 2 AC . (21) 

If we now use the factorization method from Q, with g(v) — and h(v) = X 2 v + cv~ 3 , the corresponding Abel 
equation for the EP equation (|20|) is: 

^ + h ( v ) = . ( 22 ) 

which gives 

rf = a - X 2 v 2 + cv- 2 . (23) 
We will see later that the integration constant c\ is actually an invariant of (|2"U|) . 

IV. EP EQUATION WITH ABEL-INDUCED DISSIPATION 

We build now the Ermakov-Pinney equation with Abel-type dissipation according to the equation 

V(( + g v (v)v c + h v {v) =0 (24) 

proposed by us in a previous work 0], where g v is obtained from h v using Chiellini's integrability condition for Abel's 
equation, ^ (f(Sy) = kg{u), k a real constant, which leads to 

g v = Hv . (25) 



Thus, one gets 



X 2 v 2 + cv 2 
\/kX 2 v 4 + c\v 2 — kc 
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The inverse of the Abel solution, i = jf-\k=-2, see Q, is obtained from 

- = „ " , (27) 

which gives 

C " C0 ^/7^w "w + 2c - (28) 



Let oj — y/2 Xv 2 , then (|2"8"]l becomes 



I f du 

C - Co = — =- / , == , 29 



where a = c + gw, which gives the particular solution v(C), via 



^ 2 = Cl ± V8A 2 c + ^sin(2V2A(C-Co)) (3Q) 
This equation is extremely important since the constant c\ is an invariant (Toj . 

V. GENERAL SOLUTION TO THE EP EQUATION WITH ABEL-INDUCED DISSIPATION 

We will introduce now the following Ermakov system with Abel-induced dissipation 



"CC + 9u(u)u(() + h(()u + bu 3 = 
va + 9v(v)v(C) + X 2 v + dv- 3 = . 



2„ i w„-3_n ( 31 ) 



We have the following result. 

Theorem. For the dissipative Ermakov system (|3lT) with general h(Q, and v a particular solution via ([30)l and 
given by (1261) . the general solution m(C) is obtained by 



_ r 4j ± y^T^sin (2yg(e - 9 )) d ± V8A 2 c + c\ sin (2^2A(C - Co)) 1 \ 
UtCj _ I 2^ 4A 5 J ' 

where Ibd is the following invariant 

hd = -bi K -J -d[-) + (u ( v - uv c ) 2 . (33) 

Proof: 

First, we construct the invariant Ibd, using Q33p. 

Then, we let 8 = J ^^C: an d 9 — — in (|3"3")l to obtain the following separable equation 

= d8 . (34) 



Vb + I bd 9' 2 + d6 4 



Since d = —p < 0, then (l34l) has the same format as (l28l) with p = 2A 2 , = Ci, 6 = 2c. Then, by integrating we 
obtain 



2p 

But from 9 = ^j, one gets u — 9v so that using (1531) and (j3T))) leads to 



I bd ± y/T^ + 2bpsin (2 ^(8 - 8 )) 
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We give now two related lemmas. 

Lemma 1: The invariant Im is related to the factoring junctions of the system |5I 



Let $1, and V^i be the factoring functions such that = <&\u — > u = c u e^ ld( ^, and = — > v = c v e^ ld ^. 
Then, using (|33|) we have 

I bd = 0V(* X - vKO 2 - m cosh ( j (* x - tt^C) , (36) 



/ 2 2 

where m = 2 1 + g?^- 
Lemma 2: I = C\. 

If one writes ([551) as 



using chain rule, (J35J , and by substituting O = J -^dC, leads to Ibd = c\. 
Let us choose A = c = \ , and Co = in (1501) , then 



6 = T\/2arctanh yl + cf ±Citan(C/v / 2) , (38) 

and 



9 2 = Cl ± Jc? + 1 sin(v / 2e). (39) 
By substituting both into ([57]) . with b = 1, <i = — i, and 9o = 0, then I = c\. 

VI. ABEL-DISSIPATIVE EP EQUATION WITH A 2m 1 NEGATIVE POWER NONLINEARITY 

Reid has shown in [ll[ that the general solution to 

v cc + h(C)v = q m {C)v- {2m -^ (40) 

is 

" = << + <dw*"™ )V ™ (41> 



provided that u x and u-i are two independent solutions of (|16p . and 

q m (C) = c( Ul u 2 ) m - 2 . (42) 

We are interested in a general solution to 

u a + g u (u)u c + h(C)u = q m (cy- 2m , (43) 

via the machinery of invariants. 

First, let us choose h(() = +A 2 ,0, — A 2 , then (|40|) has solutions 

w_ = (a m e mAC + c m b m e- mX<: ) 1/m , h(() = - A 2 

«o = (l + ^C m )- , MC) = o, (44) 

w+ = (a TO cos mAC + £ m 6 m sin mA() 1/m , ft(C) = A 2 , 



where 



4A 2 (ab) m {m - 1) 



(45) 
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For simplification, let us write (|44l) as 



v Q = (l + B ( m )^ 

v + = (A cos m\( + B sin toA£) 



l/r, 



(46) 



where A = a m , B = 4X i am c (m _ iy B = and by integrating 6 T = / -^d( and 9 = / -±?d( leads to 



9. 



2 

11 rn—1 . m+1 . I j^_ p ra\c\ \ 



©o = C2i ;l l(^,^; I: ^ i ;-S C , ^ 

sin 2 / m ^mA£+arctan (A/B)l cos ^mA^+arctan (A/ 



(47) 
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mA ( A cos mA(+6 sin mA( ) 



:Fi(f,| + ^;§;cos 2 (mAC + arctan (A/ B) 



Finally, the solution to (l4"5|) is u = with 9 obtained from 
Using the same constants as above, together with a = b = c = cx = l, then the solution is the following 



.(Cm) = (e+ + Jzie-^) ™ (l ± v/2sin(V29_)) 

1/2 



1/2 



u (Cm) 



1 



1 ± \/2sin(V2ec 



fe(C) = —A 2 , 

h(0 = o 



(48) 



where 



^ u+(C m) = ( cos 2* + T + I sin 2*)™(l± v / 2sin(x/29+)) 1/2 , = A 2 , 



(m 



Q _ /- JP ( J_ _2_ - m+1 . _ 1 fr. 

^ 2 llrn'm' m ' m— 1 ^ 



(49) 



e+ = 



2 sin m ^ 2|£ +arctan (rn— l)j cos (^jg^ + arctan (m— 1)1 



m cos - 



F i (s» I + h b cos2 + arctan( ? n - 1)) 



Note that for the case m = 2 which is the standard EP case the solution simplifies to 



tt-(C,2) = J{e$ + e~t) (l + V2 sin (\/2 arctan e 



< uo(C,2) = ^(C + C^ 1 ) C + V2Csin(V2arctanC) 



(50) 



u+(C,2) = J (cos C + sin C) (l + V2 sin (arctanh cos ^ c )) . 



Plots of the solutions u for different values of m and all the other parameters set to the value of unity are given in 
Figs. ([I])-®. Solutions u_ are real, while u + are periodically pure real and pure imaginary. As for uq, for m = 2, it is 
pure imaginary on the negative half-axis and pure real on the positive half-axis, while for m > 2 the solutions display 
branch cuts. Notice also the diminishing of the amplitudes with increased order of the negative-power nonlinearity. 



VII. CONCLUSION 



Ermakov-Pinney equations, eiher with standard m = 2 or higher-order m > 2 nonlinearities, have been treated as 
dissipative second-order ODEs with the dissipation induced by means of an Abel integrability criterion as discussed 
in [7|. Using the dynamic invariant of an Ermakov system built with equations of this type, the general solution of 
these equations is given in the general case of odd Ermakov nonlinearities. 
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FIG. 1: (Color online) u_(f) for the case m — 2. 
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FIG. 4: (Color online) The modulus square of it+(C) for the case m — 2. 
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FIG. 5: (Color online) The modulus square of u+(Q for the case m — 



